Forest-fire models as a bridge between different paradigms in Self-Organized 

Criticality 



Proshun Sinha-Ray and Hcnrik Jcldtoft Jensen 
Department of Mathematics, Imperial College, 180 Queens Gate, London SW7 2BZ, United Kingdom 

We turn the stochastic critical forest-fire model introduced by Drossel and Schwabl (PRL 69, 1629, 
1992) into a deterministic threshold model. This new model has many features in common with 
sandpile and earthquake models of Self-Organized Criticality. Nevertheless our deterministic forest- 
fire model exhibits in detail the same macroscopic statistical properties as the original Drossel and 
Schwabl model. We use the deterministic model and a related semi-deterministic version of the 
model to elaborate on the relation between forest-fire, sandpile and earthquake models. 



Several types of models of self-organised criticality 
(SOC) exist jUg. The original cellular automaton mod- 
els were defined by a deterministic and conservative 
updating algorithm, with thresholds (barriers to activ- 
ity), and stochastic driving |||]]. A new variation of 
models was developed by Olami, Feder and Christensen 
(OFC) H who realised that a non-conservative thresh- 
old model might remain critical if driven uniformly. The 
OFC model is completely deterministic except for a ran- 
dom initial configuration. In both types of model the 
threshold is assumed to play a crucial role as a lo- 
cal rigidity which allows for a separation of time scales 
and, equally important, produces a large number of 
metastable states. The dynamics take the system from 
one of these metastable states to another. It is believed 
that separation of time scale and metastability are essen- 
tial for the existence of scale invariance in these models. 

A seemingly very different type of model was developed 
by Drossel and Schwabl (DS) ||. No threshold appears 
explicitly in this model and the separation of time scales 
is put in by hand by tuning the rates of two stochastic 
processes which act as driving forces for the model. The 
DS forest-fire (FF) is defined on a d-dimensional square 
lattice. Empty sites are turned into "trees" with a prob- 
ability p per site in every time step. A tree can catch fire 
stochastically when hit by "lightning" , with probability 
/ each time step, or deterministically when a neighbour- 
ing site is on fire. The model is found to be critical in the 
limit p — > together with f /p —> 0. This model is a gen- 
eralization of a model first suggested by Bak, Chen and 
Tang Q which is identical to the DS model except that 
it does not contain the stochastic ignition by lightning. 
The BCT system is not critical || (in less than three 
dimensions, see A continuous variable, uniformly 

driven deterministic version [ fiof also shows regular be- 
havior for low values of p . Thus the introduction of 
the stochastic lightning mechanism appeared to be neces- 
sary, at least in two dimensions, for the model to behave 
critically. A useful review can be found in fT^j . 

In the present letter we present a transformation of 
the forest-fire model into a completely deterministic sys- 
tem. This model is an extension of the recently intro- 



duced auto-ignition forcst-firc, a simple variation on the 
DS model As in that model, we find that all macro- 
scopic statistical measures of the system are preserved. 
Specifically, we show that the three models have the same 
exponent for the probability density describing clusters 
of trees, similar probability densities of tree ages and, 
probably most unexpected, almost the same power spec- 
trum for the number of trees on the lattice as a function 
of time. It is surprising that the temporal fluctuation 
spectrum can be the same in the deterministic model as 
in the DS forest fire, since even a small stochastic ele- 
ment in an updating algorithm is known to be capable of 
altering the power spectrum in a significant way [pd[ . 

Definition of model - The SOC FF can be recast into 
an auto-ignition model. This model is identical to the DS 
model, except that the spontaneous ignition probability 
/ is replaced by an auto-ignition mechanism by which 
trees ignite automatically when their age T after incep- 
tion reaches a value T max . Choosing this value suitably 
with respect to p gives a system with exactly the same be- 
haviour and statistical properties as the DS model Jl3| . 
Thus one stochastic driving process has been removed 
and a threshold introduced, while maintaining the SOC 
state; this model also displays explicitly the relationship 
between threshold dynamics and the separation of time 
scales so necessary for the SOC state. 

The auto-ignition model can be turned into a com- 
pletely deterministic critical model by eliminating the 
stochastic growth mechanism. In the deterministic model 
(which we shall call the regen FF) each cell is given an 
integer parameter T which increases by one each time 
step. If T > 0, the cell is said to be occupied, otherwise 
it is empty (or regenerating) . The initial configuration is 
a random distribution of T- values and fires. Fires spread 
through nearest neighbours and the auto-ignition mecha- 
nism is again operative so that a tree catches fire when its 
T = T max . However, in this model when a tree catches 
fire the result is a decrement of T regen from its T-value. 
Note that when T regen < T max , a cell may still be occu- 
pied after it has been ignited. The parameters T max and 
T r egen can be thought of as having a qualitatively recip- 
rocal relationship with / and p respectively (in terms of 
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the average 'waiting time' for spontaneous ignition and 
tree regrowth), though this is less straightforward in the 
latter case because trees are not always burned down by 
fire. It is evident that T regen also sets, and allows direct 
control of, the degree of dissipation of the T-parameter 
in the system. 

Results - We now turn to a comparison between the 
statistical properties of the stochastic DS FF and the 
entirely deterministic regen model, with reference to the 
partly deterministic auto-ignition model. 

First we consider the probability density p(s) of the 
tree clusters sizes |l5| ] simulated for different parameters 
for the different models. It is well known that the corre- 
lation length in the DS model (as measured by the cut-off 
s c in p(s)) increases as the critical point is approached by 
decreasing p, / and f /p There is a corresponding in- 
crease in the power law regime for the cluster distribution 
in the auto-ignition model as p is decreased and T max is 
increased fl3| . The scaling behaviour of the cut-off s c 
is difficult to ascertain due to the limited range of data 
available, but seems to be of the form ln(iS c ) ~ pT max , 
although we cannot exclude an algebraic dependence of 
the form s c ~ (pT max ) a , with a ~ 6. Fig. |l| shows scaling 
plots for the regen model, and we see that here too the 
cut-off s c scales with increasing ratio, t — T rnax /T regen . 
We have approximately ln(s c ) ~ T max though again the 
relation may be algebraic. The conclusion is that all three 
models approach a critical state described by the same 
power law p(s) ~ s~ T with r ~ 2. 
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FIG. 1. Scaling plots of tree cluster distributions 
for regen model (L = 1000): Tmax = 200 and 
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t — 1.1. Inset: Scaling of cut-off with t increasing left to 
right (dotted) and sample distribution for DS model with 
p = 0.001, f/p = 0.01 for comparison (solid). 



the age-age correlation function: 

C(r) = (T(r + ro)T(r )}-<T(ro)} 2 



(1) 



This correlation function was never studied for the DS 
FF because the model does not consider the age T(r) 
explicitly. In Fig. || we show the behavior of the age-age 
correlation function in the regen and DS models. 
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FIG. 2. T — T correlation functions for DS (L=1000; 
p = 0.001, f/p = 0.01 - dotted) and regen (L=2000; 
T max = 200 and t = 1.0, trees and empty sites - solid and 
long dashed) models. Power law correlations are clearly ob- 
served. 

As usual it is difficult to obtain a substantial power law 
region because of finite size limitations. Nevertheless it is 
clear that C(r) does exhibit power law dependence on r 
and we find C{r) ~ r~ r > with 77 ~ 0.32, 0.21 and 0.23 for 
the regen, auto-ignition and DS models respectively. In- 
terestingly, the same correlation function for empty sites 
(which have negative T in the regen model is also a power 
law with r\ ~ 0.13. 

Let us now turn to the temporal characteristics of the 
models. In Fig. ^ we show that the probability distribu- 
tion of the ages of the trees has a very similar form for 
all three models. 

All are broad and exponential in character. Since it is 
a microscopic property, it is not surprising that there is 
some variation between the models. This variation may 
also be the reason for the different exponents in the age- 
age correlation functions mentioned above. 



One expects the power law observed in the cluster size 
distribution to be reflected in power laws for spatial cor- 
relation functions. It is particularly interesting to study 
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FIG. 3. Age profiles for DS (p = 0.001, f/p = 0.01, plot- 
ted with a T max of 4000 for clarity - dotted), auto-ignition 
(p=0.0001, T max = 20000 and T max = 24000 - dashed and 
solid) and regen (T max — 1000, t — 1.0 - long dashed) mod- 
els. 

It is remarkable that the DS FF exhibits a cut-off in 
the age distribution which is nearly as sharp as the cut- 
off in the two threshold models. This shows that the 
stochastic ignition process in the DS model, character- 
ized by the lightning probability /, can be replaced in 
surprising detail by the deterministic age threshold. 

The collective temporal behaviour is represented by 
the power spectrum of the time variation of the total 
number of trees on the lattice. In Fig. ^ these power 
spectra are shown for the DS and regen models (again, 
the power spectrum for the auto-ignition model is nearly 
identical). 
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FIG. 4. Tree density power spectra, for DS (f/p = 0.01 - 
black) and regen (T r egen/T ma x = 1.0 - grey) models 



Our most surprising result is that the deterministic 



regeneration model has nearly the same power spectrum 
as the two other models, particularly in the light of the 
differences in the age profiles above. 

The equivalence between the three models allows us to 
think of the probabilistic growth and lightning in the DS 
FF model as effectively acting as thresholds. Qualita- 
tively one can readily see that the probabilistic nature of 
the growth and the lightning can be interpreted as a kind 
of rigidity. Namely, an empty site has a rigidity against 
being turned into a tree described by l/p. A tree has 
rigidity against fire described by the fact that a tree only 
catches fire if nearest neighbor to a fire or when hit by 
lightning. 

Discussion - We now discuss the relationship between 
the regen model presented above and other SOC models. 

Our regen model is similar to the deterministic model 
introduced by Chen, Bak and Jensen [fjo|| . The crucial 
difference however, is that in the previous model the ratio 
Tregen/T m ax - which must be decreased to move closer to 
the critical point and obtain scale free behavior - is effec- 
tively held fixed at a finite value, and hence the model 
does not allow one to approach the critical state. 

The regen model has several features in common with 
the sandpile and earthquake models. It is similar to both 
sets of models in that the intrinsic dynamics is entirely 
deterministic and controlled by thresholds. The model 
is uniformly driven like the OFC earthquake model || , 
and moreover, our deterministic FF model is genuinely 
non-conservative. It is worth noting that distributing the 
increase in T randomly in a limited number of portions 
(rather than equally across all trees) each time step was 
found to destroy the criticality as the size of the por- 
tions increased. In one important respect our model is 
more similar to the BTW sandpile model than to the 
OFC model. Namely, when a site suffers relaxation (a 
tree catches fire) a fixed amount T reg en is subtracted 
from the dynamical variable of that site. The same hap- 
pens in the BTW model. In the OFC model, on the 
other hand, the dynamical variable of a relaxing site is 
reset to zero. This property has been argued to allow 
for a marginal synchronisation in the model and hence 
to be responsible for the OFC model's ability, in con- 
trast to the BTW model, to remain critical even in the 
non-conservative regime Jh| . Seen in this context the de- 
terministic FF model presented here constitutes a very 
interesting mix of features from the BTW and OFC mod- 
els. Our regen FF model is non-conservative, uniformly 
driven and though the microscopic update does not sup- 
port a marginal synchronisation, nevertheless the model 
does exhibit the same scale free behavior as the DS FF. 

This gives a direct link between the SOC behaviour of 
the BTW, OFC and DS FF models, each of which are 
commonly assumed to be representative of different and 
distinct classes of SOC models (e.g. in Further- 
more, the change in the mechanism for the renewal of 
the forest (from a probability for growth to a time for 
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regeneration) and the resultant sandpile-like picture al- 
lows the identification of p with a dissipation parameter 
(in terms of the subtraction of T regen on ignition) rather 
than as a driving parameter. This is quite contrary to 
the normally held and most obvious view - for the DS FF 
- that p is the driving parameter (creating trees in the 
system), and that if anything / controls the dissipation 
(the complete combustion of trees into empty sites). If 
this is so, we can speculate that it may be possible to re- 
late the physical limits for critical behaviour in the BTW 
sandpile: 

h, h/e -> 

(where h is the driving rate and e the dissipation) and, re- 
calling the qualitatively reciprocal relationships between 
f,P,T max and T regen noted earlier, the DS and regen 
forest-fire models: 

/? f /p y 0, and ^/T max: T re g en /T max ► 

The main difference between the deterministic FF 
model and the sandpile and earthquake models is that 
the dynamical variable T is not transported to neighbor- 
ing sites when a site relaxes and that the threshold exists 
only for the initiation and not propagation of avalanches. 
This difference can be summarized as the FF model being 
a model of two coupled fields, fires and trees, whereas the 
sandpile and earthquake models contain one self-coupled 
field, the energy of a site. 

Another difference consists in that the thresholds of 
the deterministic FF model must be tuned (to infinity) 
for the model to approach the critical regime. The reason 
for this is that the thresholds relate directly to the rate 
of driving in the model. The sandpile and earthquake 
models are different in that the SOC limit of slow driving 
can be reached without a tuning of the thresholds. 

Finally, we note that the regen model is critical with 
periodic boundary conditions (in contrast to the BTW 
and OFC models) and without external driving (unlike 
the DS model), and is therefore the only system which 
can be said to be completely self-contained. 

Conclusion - We have demonstrated that the stochas- 
tic Drossel-Schwabl forest-fire model can be turned into 
a deterministic threshold model without changing any 
of the collective statistical measures of the system in a 
significant way. The model illuminates greatly the rela- 
tionship between different types of SOC models. 
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